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ABSTRACT 

We have newly identified a substantial number of type 1 active galactic nuclei (AGN) featuring weak 
broad-line regions (BLRs) at z < 0.2 from detailed analysis of galaxy spectra in the Sloan Digital Sky 
Survey Data Release 7. These objects predominantly show a stellar continuum but also a broad Ha 
emission line, indicating the presence of a low-luminosity AGN oriented so that we are viewing the 
central engine directly without significant obscuration. These accreting black holes have previously 
eluded detection due to their weak nature. The new BLR AGNs we found increased the number of 
known type 1 AGNs by 49%. Some of these new BLR AGNs were detected at the Chandra X-ray 
Observatory, and their X-ray properties confirm that they are indeed type 1 AGN. Based on our new 
and more complete catalogue of type 1 AGNs, we derived the type 1 fraction of AGNs as a function of 
[O III] A5007 emission luminosity and explored the possible dilution effect on the obscured AGN due 
to star-formation. The new type 1 AGN fraction shows much more complex behavior with respect to 
black hole mass and bolometric luminosity than suggested by the existing receding torus model. The 
type 1 AGN fraction is sensitive to both of these factors, and there seems to be a sweet spot (ridge) 
in the diagram of black hole mass and bolometric luminosity. Eurthermore, we present a hint that the 
Eddington ratio plays a role in determining the opening angles. 

Subject headings: galaxies: active — galaxies: Seyfert — galaxies: nuclei — quasars: general — 
galaxies: statistics — methods: data analysis 


1. INTRODUCTION 

According to the unification scheme of active galactic 
nuclei (AGNs, Antonucci 1993; Urry & Padovani 1995), 
which is widely accepted in current astrophysics, differ¬ 
ent viewing angles explain the variety of phenomenologi¬ 
cal subclasses of AGNs. For example, the optical spectra 
of type 1 AGNs present broad permitted emission lines 
as well as narrow emission lines exposing their central 
nuclear region, whereas type 2 AGNs show only narrow 
emission lines. According to the simplest AGN unifica¬ 
tion model and its strictest interpretation, the feasibil¬ 
ity of observing broad permitted emission lines depends 
solely on the angle between the observer’s line of sight 
and the axis of an obscuring medium, often called the 
“dust torus”. 

With the aid of mega-scale surveys such as the Sloan 
Digital Sky Survey (SDSS; York et al. 2000) that provide 
homogenous photometric and spectroscopic data, AGNs 
have become one of the most actively studied topics in 
astrophysics. In particular, studies on obscured AGNs 
showing narrow emission lines have shed light on the con¬ 
nection between the host galaxies and the central AGN 
(Kauffmann et al. 2003; Salim et al. 2007; Schawinski 
et al. 2007; Westoby et al. 2007; Schawinski et al. 2009; 
Haggard et al. 2010; Schawinski et al. 2010). 

Unobscured (type 1) AGNs reveal their bright nuclear 
light directly, and it has been an important task and 
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subject of interest to decompose the observed light into 
a central nuclear component and a diffuse galaxy compo¬ 
nent. For example, the high-spatial-resolution Advanced 
Gamera for Surveys on the Hubble Space Telescope, has 
been extensively used for this purpose (Sanchez et al. 
2004; Ballo et al. 2007; Alonso-Herrero et al. 2008; Sim¬ 
mons & Urry 2008; Schawinski et al. 2011; Simmons 
et al. 2013) in combination with two-dimensional surface 
brightness fitting techniques (Peng et al. 2002). 

The fraction of type 1 AGN (among all AGNs) can 
act as an effective test of any model for AGN classifica¬ 
tions and can be used to extract and constrain critical 
information, such as the opening angle and source lumi¬ 
nosity, in the assumed framework. Based on a sample 
of high-luminosity and radio-selected broad-line AGNs, 
Willott et al. (2000) found the type 1 AGN fraction of 
40% and a typical value of the “half-opening” angle of 
about 53°. The precise value of the type 1 AGN fraction 
has been reported to be substantially lower than this 
finding, when low-redshift surveys of lower optical lumi¬ 
nosity cuts are used instead (Osterbrock & Shaw 1988; 
Huchra & Burg 1992; Maiolino & Rieke 1995; Maia et al. 
2003). This difference seems to be sensible, as more re¬ 
cent, multiwavelength studies found a strong, positive 
luminosity dependence of the type 1 AGN fraction (Stef¬ 
fen et al. 2003; Hao et al. 2005; Simpson 2005; Hasinger 
2008; Burlon et al. 2011; Assef et al. 2013). 

AGNs exhibit characteristic spectral features across 
the electromagnetic spectrum: e.g., [Oiii] A5007 optical 
emission lines, mid-infrared and hard X-ray luminosity as 
isotropic quantities tracing the source luminosity, which 
are naturally employed by many type 1 AGN fraction 
studies. In particular, Simpson (2005) found that the 
type 1 AGN fraction gradually increases with L[Oiii] 
based on SDSS Data Release 2 (Abazajian et al. 2004) 
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Fig. 1. — Flux ratio (F6,533/F6,47o) distribution versus signal-to-statistical noise (S/sN). Demarcation lines in the central panel are 
shown as red dashed lines. The three lines indicate 84.1%, 97.7%, and 99.8% distributions. Example spectra represented by the numbered 
markings in the central panel are shown in the surrounding panels, with two small sub panels for each object. Surrounding panels marked 
with ‘a’ on the top left side show the overall shape of the spectrum, whereas the panels marked ‘b’ show the Ha region. Each object is 
labeled in sub panel ‘a’ by SDSS ObjID. Black solid lines indicate the observed spectrum. The two narrow bands that define the flux ratio 
are shaded gray. Four spectra (numbered from 1 to 4) are chosen as type 1 AGN candidates out of six examples by the demarcation lines. 
Flux has the unit [10“^^ergs~^cm“^A~^]. 


and hard X-ray data in the literature (Ueda et al. 2003; 
Grimes et al. 2004; Hasinger 2004). 

A theoretical framework has been claimed to explain 
such a luminosity dependency of the type 1 AGN frac¬ 
tion. For example, the receding torus model (Lawrence 
1991; Falcke et al. 1995; Hill et al. 1996; Simpson 2005) 
imposes a dust sublimation radius of the torus, which is 
directly affected by the nuclear luminosity. As the lumi¬ 
nous central engine pushes the dust torus away, thereby 
increasing the torus’ inner radius, the opening angle must 
increase if the height of the torus is fixed. This naturally 
explains the observed luminosity dependence of the type 
1 AGN fraction. 

A precise determination of the type 1 AGN fraction is 
obviously a key test of a model’s validity. Since the AGN 
unified model explains the observability of both types of 
AGN as the consequences of obscuration induced by the 
geometry and the structure of AGN, the frequency of ob¬ 
served type 1 AGN indirectly constrains the model. For 
instance, higher type 1 AGN fraction implies wider open¬ 


ing angle, smaller scale height of the torus, and higher 
nuclear luminosity. Herein we report on our efforts to 
discover new type 1 AGNs from the SDSS database to 
provide such a test. We describe how we found a sub¬ 
stantial number of new type 1 AGNs and derived a more 
up-to-date type 1 AGN fraction. The use of the new type 
1 AGN fraction places a harsh constraint on the model, 
and if the scheme used herein proves convincing, it ap¬ 
pears that an important modification to the model will 
be necessary. We release a new catalogue of type 1 AGNs 
in the nearby universe by means of this paper. 

The paper is organized as follows. In Section 2, we 
describe the selection technique whereby we constructed 
the new type 1 AGN database, and describe a simulation 
that quantified the detection efficiency of our method. 
We demonstrate our conservative scheme for type 1 se¬ 
lection and determine the final type 1 AGN sources, tak¬ 
ing the simulation into account. In Section 3, we present 
the new type 1 AGN fraction and discuss how contami¬ 
nation caused by star formation activities may affect it. 
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Finally, we discuss our findings and summarize our re¬ 
sults in Section 4. 

We assume a cosmology with h = 0.70, Um = 0.30, 
and Ua = 0.70 throughout this work. 

2. TYPE I AGN SAMPLE SELEGTION 

2.1. Selection based on broad Ha features 

The presence of unidentified type 1 AGNs in the nearby 
universe has been reported by Hao et al. (2005), Greene 
& Ho (2007), and Oh et al. (2011) [a.k.a. the OSSY 
catalogue]. In particular. Oh et al. (2011) found that the 
presence of broad Balmer emission lines were responsible 
for a very poor fit by model to both the emission lines and 
the stellar continuum during their construction of the 
OSSY catalogue. The authors used a flux ratio near the 
Hof emission line to systematically search for BLR AGNs 
from the SDSS DR7 galaxy database. In this work, we 
adopted their basic scheme but developed it further to 
minimize dubious detections. 

The OSSY catalogue provides new spectral line 
strength measurements for the SDSS DR7 galaxies at 
z < 0.2 (N = 664,187). By combining stellar kinemat¬ 
ics measurements (pPXF; Gappellari & Emsellem 2004) 
and spectral line fitting algorithm (gandalf; Sarzi et al. 
2006), the authors released improved stellar absorption 
and emission-line strengths, stellar velocity dispersion, 
and quality assessing parameters. We start our analysis 
based on the OSSY catalogue. 

To identify the spectra showing broad Ha features®, 
we computed a ratio between the mean fluxes at two 
specific wavelength bands. The first narrow band (6,460 
— 6,480 A) is located near Ha but far enough from it to 
be free from the Ha broad emission, mainly to secure the 
pseudo-continuum for Ha. The second (6,523 — 6,543 A) 
was placed close to [N ii] A6548 to highlight the broad Ha 
feature. Then, we computed a flux ratio (Eg 533 /E 6 , 47 o) 
based on the two bands. Fig. 1 shows the flux ratio versus 
the signal-to-statistical noise ratio (S/sN) of the chosen 
specific wavelength bands for all galaxies listed in the 
OSSY catalogue (N = 664,187; z < 0.2). We used this 
diagram to first detect Ha BLR candidates. Examples of 
spectral energy distributions (SEDs) for various values of 
flux ratios are illustrated in Fig. 1 (surrounding panels), 
including the overall spectra and the Ha regions. As 
the flux ratio increases, broad Ha features become more 
prominent (see objects marked 1, 3, and 5 in order of 
decreasing broad Ha strength), and the same trend is 
visible in the higher S/sN regime (objects 2, 4, and 6). 
It seems clear that the flux ratio is effective for detecting 
Ha BLR candidates. 

The Icr demarcation line shown in Fig. 1 was applied 
to first select type 1 AGN candidates. A total of 17% 
of galaxies have a flux ratio above this cut. We used 
gandalf (Gas and Absorption Line Fitting, Sarzi et al. 
2006) to fit the spectrum based on both nebular emis¬ 
sions and stellar absorptions. Our approach for analyzing 
the SDSS spectra can be summarized as a two-step pro¬ 
cess. First, we masked nebular emission lines and derived 
the stellar kinematics by directly matching the observed 
spectrum with stellar templates by using the pPXF IDL 

® broad Hq is also a signature of supernovae but it is known as 
rare (Filippenko 1997). 


TABLE 1 

Simulation designs 


components 

range 

increment 

No. of bins 

S/sN^ 

5 - 55, > 55 

5 

11 

FWHM 

50- 15,000 [km s”*] 

50 

300 

Peak*’ 

0.01 X — 0.09 X Ahck 

0.01 

9 

0.1 X Alia ~ 1-0 X Alia 

0.1 

10 


Note. — In total, 627,000 mock galaxies were generated. 

^ Average S/sN from 4 continuum bands (4,500-4,700, 5,400-5,600, 6,000- 
6,200, 6,800-7,000 A) 

^ Peak of the Gaussian 
^ Amplitude of the narrow Ha 

program (Gappellari & Emsellem 2004). Eor the “stel¬ 
lar” templates, we used the Bruzual & Gharlot (2003) 
stellar population synthesis models and the MILES stel¬ 
lar library (Sanchez-Blazquez et al. 2006). After measur¬ 
ing the kinematic broadening, we lifted the masking on 
the nebular emission lines and simultaneously fit stellar 
continuum and emission lines by using the Levenberg- 
Marquardt minimization (Markwardt 2009)^. 

The standard pPXF and gandalf procedures generally 
work well on normal galaxies, but a more elaborate con¬ 
sideration is necessary for type I AGNs. For example 
only for type I AGNs, we adopted a wide masking width 
of 12,000 km s“^ for broad Balmer lines (Ha, Hy, 
and H(5) and 1,200 km s“^ for narrow emission lines, 
prior to measuring stellar kinematics. Besides, we im¬ 
posed an initial guess of a large width, 1,000 km s“^, for 
broad Balmer component fitting. We implemented sin¬ 
gle broad Gaussian component for each Balmer line. As 
demonstrated in the middle panels of Fig. 2, these pre¬ 
scriptions were effective for decomposing the broad (Ha) 
and narrow (Ha, [N ii] A6548 and [N ii] A6584) compo¬ 
nents around the broad Ha region. 


Following these routines, we measured the full-width at 
half-maximum (FWHM), flux, and amplitude over noise 
(A/N) ratio® of all lines listed in the OSSY catalogue. 
Our selection criteria for type 1 AGN candidates are as 
follows. 

• 0.00 < z < 0.20 

• FWHM of broad Ha > 800 km s“^ 

• A/N of broad Ha > 3 

The redshift limit is constrained by the wavelength of 
the red end of the pseudo-continuum around Ha. We 
chose a cut for FWHM(Ha) that was slightly lower than 
the one more widely used, 1,000 km s“^ (Vanden Berk 
et al. 2006; Schneider et al. 2010; Stern & Laor 2012), 
hoping to detect even weaker BLR AGNs. As a result, 
we found 9,671 (8.6%) type 1 AGN candidates among 
the 111,824 sources that satisfied the Icr cut shown in 

^ pPXF and gandalf use the MPFIT IDL routine provided by 
Craig Markwardt at http://cow.physics.wise.edu/'craigm/idl/ 

® A/N is defined by the ratio between Gaussian amplitude of 
the emission line and dispersion of continuum residual of the given 
spectrum 
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Fig. 2. — Example spectra with SDSS gri composite images for a newly identified type 1 AGN (top) and a previously known type 1 
AGN from the SDSS pipeline (bottom). Left and middle panels present the overall spectra and their Hct regions, respectively. Black lines 
indicate the observed spectra and red dashed lines include stellar continuum fits. Note that the observed spectrum and the fit in the top 
panel are shown with offset in the ordinate for clarity (+50 in flux unit). Detected emission lines are shown in the left panels as blue lines. 
In particular, broad Ha components are shown in the middle panels as well as other narrow emission lines. Dots in the middle panels show 
the residuals. 
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Fig. 3. — Flux ratio distribution with respect to S/sN (left) and FWHM (right) for the simulated galaxies. In the left panel, color grids 
and the corresponding color bar represent the mean FWHM. In the right panel, FWHMs measured for the actively star-forming galaxies 
are indicated by a blue contour (95%) while the values for the mock spectra are indicated by achromatic filled contours. 
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TABLE 2 

Detection efficiency for each a cut 


cr CUt^ 

No. of 

mock galaxies 

No. of 
mock BLRs 

completeness^ 

(%) 

purity 

(%) 

> Icr 

547,495 

543,052 

91.49 

99.19 

> 2cr 

515,939 

514,562 

86.69 

99.73 

> 3(7 

487,943 

487,512 

82.13 

99.91 


^ cr cut marked in the central panel of Fig. 1 

^ Number of all mock broad-line AGNs (FWHM > 800 km s“^) — 
593,560 


TABLE 3 

Detection efficiency for each a cut after applying 
AREAL RATIO 


(7 cut 

No. of 

mock galaxies 

No. of 
mock BLRs 

completeness 

(%) 

purity 

(%) 

> Icr 

547,495 

484,791 

81.68 

88.55 

> 2cr 

515,939 

476,341 

80.25 

92.33 

> 3(7 

487,943 

463,475 

78.08 

94.99 


the central panel of Fig. 1. Of these 9,671, only 120 
were of FWHM (broad Ha) < 1,000 km s“^, and hence 
the choice of the 800 km cutoff instead of the more 
typical 1,000 km s“^ cutoff influenced the results and 
conclusions little. 


2.2. Test of the selection scheme 

The selection process described above relies on the spe¬ 
cific flux ratio as a proxy to the broad Ha feature. To 
verify our approach and measure the detection efficiency, 
that is, its completeness and purity, we conducted a sim¬ 
ulation by generating mock spectra, each with an artifi¬ 
cially broadened Ha component. 

The basis of the mock spectra was 110 star-forming 
galaxies with narrow emission lines: 10 galaxies ran¬ 
domly chosen from each of 11 S/sN bins. To secure a 
wide range of the amplitude of Ha, we used the spectra of 
actively star-forming galaxies of SFR > 0.5 M 0 yr“^ that 
were identified by the HPT diagnostics diagram (Bald¬ 
win et al. 1981) and the Kennicutt law (Kennicutt 1998). 
We then add broad line components by varying two pa¬ 
rameters: Gaussian peak and FWHM, while the center 
of the Gaussian was fixed at 6562.8 A. The peak of the 
Gaussian was fractionally scaled by the amplitude of the 
narrow Ha emission line (Ahq) in 19 steps, namely, by 
1 % increments from 1 % to 10 % and a 10 % increment 
from 10% to 100% in Ana- The FWHM of the broad 
Ha line was varied between 50 and 15,000 km s“^ in 50 
km s“^ increments. For reference, the typical range of 
FWHM(Ha) of star-forming galaxies is 50—350 km s“^ 
within the 2cr level. In total, 627,000 (110 x 19 x 300) 
mock spectra were generated, as summarized in Tab. 1. 


Fig. 3 shows flux ratios versus S/sN (left panel) and 
FWHM (right panel) for the simulated galaxies. In 
the left panel, the galaxies with low FWHMs (dark red 


color) were generally found below the Icr demarcation 
line and/or had low S/sN. As the color grids show, our 
1(7 cut for recovering BLRs with a large value of FWHM 
(blue color) seems effective, although it still suffers from 
contaminants such as objects having noise fluctuations or 
having a broad Ha component that is too shallow. We 
will discuss how these contaminants can be removed in 
the next section. The right panel of Fig. 3 shows that the 
highly star-forming galaxies having only narrow emission 
lines obviously occupy the region of low flux ratio (< 1 . 1 ) 
and low FWHM (< 500 km s“^). Below the FWHM of 
500 km s“^, the mock galaxies have almost the same flux 
ratios as narrow emission line galaxies, which means that 
it is not easy to distinguish between weak type 1 AGNs 
and star-forming galaxies. On the other hand, the flux 
ratio of the mock galaxies begins to spread out greatly 
for FWHM > 1, 000 km s-k 

Detection efficiency (completeness and purity) can 
be measured by assuming that a BLR AGN has 
FWHM(broad Ha) > 800 km s“^. Completeness is de¬ 
fined as the fraction of mock BLR AGNs that are recov¬ 
ered by the procedure (i.e., the a cut). Purity describes 
the fraction of BLR AGNs among all galaxies selected by 
a cr cut. Tab. 2 lists the completeness and purity of BLR 
AGNs for three choices of tr cut. As higher cr cuts are 
applied, completeness decreases but purity increases. 

Fig. 4 illustrates that the completeness (color key) of 
our selection method based on the flux ratio of the broad 
Ha component is mainly determined by the broad Ha 
luminosity (ordinate). We used this completeness distri¬ 
bution to find BLR AGN candidates most economically. 
We first located the locus of 90% completeness pixels 
in the diagram. As most of the pixels above this locus 
satisfy the completeness of 90%, the total completeness 
when using this 90% cut is in practice much greater than 
90%. Roughly 21% (2,052 out of 9,671) of the BLR AGN 
candidates (by the Itr criterion) were below the 90% de¬ 
marcation line and thus were missed by this scheme. 

2.3. Selection based on the broadening around 
[NH]A6584 

We have so far described our selection criterion based 
on the broad Ha feature, but there still are complexities. 
For instance, an extremely shallow broad Ha component 
can be technically classified as type 1 AGN as long as it 
satisfies our conditions on A/N and FWHM but not with¬ 
out doubt. For another example, the credibility of type 1 
classification for border-line objects of FWHM(Ha) just 
over 800 km s“^ can be doubted, because it is difficult to 
measure the width of such a marginally-broad Ha com¬ 
ponent accurately, especially in the midst of strong neigh¬ 
boring narrow emission lines ([N ii] A6584). 

Having carried out a sophisticated spectral fit for each 
galaxy spectrum and having determined decomposed 
broad Ha properties as a result, we were able to per¬ 
form a more stringent type 1 classification than carried 
out so far. We compared the area of the broad Ha com¬ 
ponent in question with the noise level of the observed 
continuum. The area of broad Ha beyond [N ii] A6584 to 
longer wavelengths can be used as an alternative measure 
to FWHM for the broad Ha feature. Meanwhile, the area 
defined by the dispersion (measurement uncertainty) of 
the continuum can be used as a control parameter. Fig. 5 
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Fig. 4. — Luminosity versus FWHM for broad Ho components. Completeness derived from a simulation is indicated by the underlying 
color grid. Labeled black dotted lines represent Mbh in units of solar mass following the formulae developed by Greene & Ho (2005). The 
red dashed curve indicates the demarcation used to choose type 1 AGN sources based on a 90% completeness criterion. The type 1 AGNs 
(N=7,619) selected by the 90% completeness cut are shown as blue-filled dots; excluded points (N=2,052) below the curve are shown as 
gray-filled dots. 


provides an illustration with example spectra. The right 
panels demonstrate how we assess the broad Her compo¬ 
nent (gray shaded area) that surpasses the level of the 
continuum dispersion (hatched area), which is marked as 
o'cont.- As an additional type 1 AGN criterion, we take 
the areal ratio of the broad Ha over the measurement 
uncertainty around [N ii] A6584 greater than 2. The ap¬ 
plication of this new cut removed 3,940 controversial type 
1 AGN candidates (Fig. 6). Tab. 3 lists the completeness 
and purity for each a cut after implementing areal ratio 
analysis. As is the case with any other criterion, some of 
the objects removed in this step are likely to be true type 
1 AGN. Yet, we felt it necessary to apply this additional 
criterion to construct a more robust list of type 1 AGN 
candidates. 

By using the completeness (> 90%) and areal ratio 
(> 2) arguments, we identified 5,553 type 1 AGN can¬ 
didates. It should be noted that the majority (91.3%) 
of the type 1 AGN candidates below the 90% complete¬ 
ness cut in Fig. 4 were eliminated through the areal ex¬ 
cess test, confirming the effectiveness of the completeness 
test as a filtering scheme. The flowchart in Fig. 7 shows 
the number of objects that survived each selection. The 
broad nature of the Ha emission line of type 1 AGNs 
eluded our initial attempt to measure their strengths. 


Perhaps for the same reason, many (38%, 700 out of 
1,835) of our newly found type 1 AGNs are not included 
in the MPA-JHU catalogue®. Hereinafter, our analysis is 
based on these 5,553 sources. 

We list the measured properties of these sources as 
well as the estimated quantities with designation flag in 
Tab. 9. We used the formula developed by Greene & Ho 
(2005) that depends on the line width and luminosity of 
the broad Ha to estimate black hole mass. We assumed 
bolometric luminosity by adopting the relationship de¬ 
veloped by Heckman et al. (2004) (Lboi ~ 3,500L[0ill]). 
The cross-matched flag presents that 456 type 1 AGNs 
(8.2%) and 2,954 (53.2%) sources in our catalogue are 
found in Hao et al. (2005) and Greene & Ho (2007), re¬ 
spectively. 


2.4. SDSS SpecClass=3 objects 

SDSS DR7 provides a spectroscopic classification of 
galaxies: SpecClass. Broad-line AGNs, especially for 
low redshift, are flagged SpecClass=2> when the observed 
spectrum satisfies the following criteria: the FWHM of 

® http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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Fig. 5.— Example spectra illustrating the conservative selection criterion based on areal excess. The overall spectra are shown in the 
left panels, including SDSS ObjID labels at the top left. Black and red lines respectively indicate the observed spectrum and our fit. The 
regions that are used to measure the averaged continuum dispersion are shaded gray. The middle panels present the detected [N ii] A6548, 
narrow Ha, broad Ha, and [Nil] A6584 emission lines (blue Gaussians), and their residuals (black dots). The right panels enlarge the 
middle panels, comparing the measured broad Ha exceeding [N ii] A6584 (gray shades) to the level of the continuum dispersion (hatched 
regions with black horizontal lines). The areal excess, that is, the ratio between the gray shaded and hatched regions, is marked, as well as 
the level of the continuum dispersion, <Jcont.- 
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N=664,187 


N=lll,824 


N=9,671 


N=7,619 


N=5,553 



(see Fig. 1) 


(see Fig. 4) 


(see Fig. 5 & 6) 


Fig. 6.— Histograms of area ratio for all type 1 AGN candidates 
(black), newly found type 1 AGNs (blue), and SDSS SpecClass 3 
sources (red). The new selection criterion, area ratio > 2, is marked 
using a vertical dashed line; it removes 2,066 dubious candidates 
(gray shaded area). 


TABLE 4 

Detection of SpecClass 3 sources 



N 

% 

total 

4125 

100.00 

detected®’ 



after flux ratio cut 

4079 

98.89 

after A/N^hq & FWHMbHa cut 

3969 

96.22 

after completeness cut (>90%)^ 

3921 

95.05 

ater areal excess cut (>100%)^ 

3718‘i 

90.13 

missed 

407 

9.87 

Description on the sources that failed flux ratio cut 


only narrow emission lines 

19 

41.30 

absence of spectral information 

11 

23.91 

maybe broad-line AGNs 

10 

21.74 

passing A/NbHa & FWHMbHa cut 

(10) 

(21.74) 

passing completeness cut (>90%) 

(8) 

(17.39) 

passing areal excess cut (>100%) 

(5) 

(10.87) 

asymmetric bump on Ho 

3 

6.52 

non-emission lines 

2 

4.35 

unknown spectrum 

1 

2.17 


^ Number of all detected SpecClass 3 sources chosen from the flux ratio 
analysis at 2 < 0.2 
^ see Section 2.3 and Fig. 4 
^ see Section 2.4 and Fig. 5 

^ Number of SpecClass 3 sources in the flnal catalog 


any emission line is greater than 1,000 km s“^, the height 
of the emission line is at least 3 times higher than its noise 
(equivalent to our condition of A/N > 3), and the equiv¬ 
alent width of the line is larger than 10 A. In total, there 
are 4,125 SpecClass 3 sources at z < 0.2. Galaxies not 
showing broad-emission lines are classified as SpecClass 
2 . 

We tried applying our type 1 AGN selection technique 


Fig. 7.— Flowchart of the type 1 AGN selection. 


588017567093293191 

_ 




. 1 . . 



4000 4600 5000 5500 _ 6000 6500 6400 6500 6600 6700 

wavelength[/l] wavelength[5l] 



wavelength[A] wavelenglh[A] 



4000 4500 5000 5500 ^ 6000 6500 6400 6500 6600 6700 

wavelength[/l] wavelength[Jl] 



wavelength[A] 


wavelenglh[A] 


Fig. 8.— Spectra of example SpecClass 3 sources below the la 
demarcation line, for a narrow-emission-line galaxy (top), an in¬ 
complete spectrum near Ho (upper), an asymmetric bump feature 
(lower), and an unknown source (bottom). Right panels show the 
Ho; region; the two narrow bands that define the flux ratio are 
shaded gray. 
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Fig. 9. — Histograms of FWHMbHa> EWbHcti logt/bHai and logL^ol for newly found type 1 AGNs (hatched area, see Section 2.3), 
SpecClass 3 sources (solid lines), and all type 1 AGN candidates (dotted) that were chosen by measuring the flux ratio F 6 , 533 /F 6,470 (see 
Section 2.1). 


on the SpecClass 3 objects; Tab. 4 summarizes our re¬ 
sults. To begin with, 4,079 out of 4,125 objects passed 
our flux ratio test. Then, 3,969 of the remaining objects 
passed the A/N and FWHM cuts, and 3,718 (90.13%) 
passed the final areal excess test. Overall, our technique 
rejected roughly 10% (407 out of 4,125) of the objects. 

We describe further why some of the SDSS SpecClass 
3 objects are not classified the same way in our analysis 
based on a subsample of the objects that were excluded 
by the first and most important step of our scheme, the 
flux ratio test. By performing additional spectral fitting 
to the 46 (4,125 — 4,079) sources that failed to pass the 
Icr cut from the flux ratio distribution, we found that 11 
objects have empty spectral bins around Ha (see Fig. 8, 
second row). Hence, our technique could not be feasibly 
applied to these objects. In the majority of the rest 35 
missed objects (25 out of 35), we could not detect any no¬ 
table broad features; the top row of Fig. 8 shows such an 
example. Of the 46 missed objects, 10 passed our A/N 
and FWHM tests of Ha and thus possibly have broad 
Ha components. However, 2 of the 10 failed to pass our 
90% completeness test and a further 3 failed to pass the 
areal excess cut test; overall, only 5 out of 10 passed all of 


our test. Similarly, our technique could not be effectively 
applied on 3 objects because their spectrum around Ha 
appeared asymmetric and uneven (see Fig. 8, third row). 
Finally, one object had a spectrum that did not appear 
to be an ordinary galaxy spectrum, making it difficult 
to fit this spectrum with a combination of template stel¬ 
lar (population) spectra (Fig. 8, bottom row). Tab. 4 
summarizes the statistics of the SpecClass 3 sources. 

Fig. 9 presents the distributions of the widths and lu¬ 
minosities of the broad Ha component for the newly 
found type 1 AGNs (hatched area) and SpecClass 3 
sources (solid line). The new sources show broad Ha 
widths (FWHMbHa) comparable to those of the Spec- 
Class 3 sources. On the other hand, the characteristic 
luminosity (i.e., the peak of the distribution) and the 
equivalent width of their broad Ha components seems 
lower, by a factor of a few, than that of the SpecClass 
3 objects, where bolometric luminosity is derived from 
Lboi ~ 3,500T[Oiii] (Heckman et al. 2004). Hence, the 
newly found objects are mostly low-luminosity type 1 
AGNs. 
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6450 6500 6550 6600 6650 6700 

wavelength[X] 




Fig. 10.— Example Chandra X-ray data (obsid = 2779) and model fit, with the SDSS gri composite image and our spectral fit. Top left: 
SDSS gri composite image. Top right: observed spectrum (black), our fit (red), detected emission lines (blue) including broad components, 
and residuals (black dots). Note that the detected emission lines are arbitrarily shifted towards the y-axis for clarity. Bottom left: Chandra 
X-ray image, including source (green circle) and background (red annulus) markers. Bottom right: extracted data with the applied fit 
versus energy (top); residuals (bottom). 


TABLE 5 

X-RAY FITTING STRATEGIES 


count 

fitting prescription 

N-ti [cm-^]^ 

rb 

< 50 

lO^'^ fixed'^ 

1.9 fixed 

50< count <100 

]^g20-24 variable 

1.9 fixed 

> 100 

]^g2u-24 variable 

1.5-2.5 variable 


^ Hydrogen column density 
^ X-ray photon index 
^ The value is fixed during spectral fit 


2.5. Chandra X-ray data 

A strong X-ray flux is one of the prominent and char¬ 
acteristic features of an AGN, because an accretion disk 
surrounding a central black hole produces X-ray emis¬ 
sions (Fabian et al. 1989). Relativistic electrons from the 
corona in the vicinity of the accretion disk scatter the 
lower-energy photons radiated from the accretion disk. 
As a result, the lower-energy photons receive more energy 
from the electrons (inverse Compton scattering), and this 


process produces an X-ray spectrum with a power-law 
continuum (Haardt & Maraschi 1991; Zdziarski et al. 
2000; Kawaguchi et al. 2001). Thus, strong X-ray emis¬ 
sion is considered to be evidence of the presence of an 
AGN. 

We utilized the Chandra Source Catalog (Evans et al. 
2010) and archival Chandra X-ray data to identify X-ray 
sources among our AGN candidates. We downloaded the 
matched X-ray images from the Chandra archive and ap¬ 
plied a 10 arcmin matching radius. Through visual in¬ 
spection of X-ray images of their optical counterpart, we 
found 84 X-ray sources among the type 1 AGN samples. 
The X-ray recovery rate (1.5%) was very low mainly be¬ 
cause X-ray sky coverage was much narrower (by a factor 
of 25) and depth was much shallower compared to the 
SDSS optical survey. 

We used the Chandra Interactive Analysis of Obser¬ 
vations (CIAO) v4.6^*’ and XSPEC vl2.8.1g software^^ 
(Arnaud 1996) to reprocess the Chandra X-ray data and 

http://cxc.harvard.edu/ciao/ 

^^ http://heasarc.nasa.gov/xanadu/xspec/ 
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TABLE 6 

X-RAY PROPERTIES OF TYPE 1 AGNS 


SDSS ObjID 

obsid'^ 

RA^ 

[deg] 

DEC‘> 

[deg] 

z 

counts 

count rate 
[count s~^] 

logi'2-lOkev'^ 
[erg s-l] 

note*^ 

587738372207476987 

3784 

121.94200 

21.23800 

0.14215 

230 

0.012 

42 yn+O OS 
<o_Q 06 

3 

587735665842520176 

2779 

211.12100 

54.39800 

0.08333 

318 

0.022 

^z.oD_o 06 

N 

587725551741370430 

2033 

143.96500 

61.35300 

0.03933 

332 

0.007 

41.36t°;°l 

N 

587741602568798223 

4695 

184.53300 

28.17600 

0.17773 

449 

0.045 

43.30^^5 

N 

587725551197683729 

7146 

124.83400 

50.00700 

0.13969 

171 

0.022 

'-0.09 

3 

587729774220738649 

9557 

210.21899 

-1.75300 

0.14878 

1662 

0.034 

42.92t°°l 

N 

587731521745518723 

827 

140.28600 

45.64900 

0.17448 

4617 

0.246 

44.16t°;°? 

3p 

587729653959033079 

887 

249.20200 

41.02900 

0.04737 

287 

0.004 

4l.28l[;;°5 

N 

587735696453009410 

1623 

220.76199 

52.02700 

0.14121 

2851 

0.190 

43.86t“;°= 

3p 


Note. — (This table is available in its entirety in a machine-readable form in the online journal. A portion is shown 
here for guidance regarding its form and content.) 

^ Chandra observation ID 
^ Chandra J2000 coordinates 

^ Errors are estimated with the 90% confidence region for a single interesting parameter. 

^ Discriminating flag for the newly identified type 1 AGNs (‘N’) and the previously confirmed sources (‘3’). Notation 
’p’ indicates the photon pileup fraction greater than 10%. 

perform a spectral fitting analysis. We manually aligned 
the coordinates to the actual image center and deter- 
mined the source aperture, Ag, (green circle in Fig. 10) 
by reading off the value from the pipeline relation be¬ 
tween the off-axis angle and the aperture size suggested. 

The applied aperture (Ag = 11.2kpc, cta^ = lO.Okpc) 
generally covers the whole galaxy (PetroROOr = 12.2kpc, 
cpetroRgOr = 6.2kpc, where PetroROOr is the radius con¬ 
taining 90% of Petrosian flux in r-band) while some 
nearby sources represent central part. Then, we used 
specextract with psfcorr in CIAO which incorporates 
an energy-dependent aperture correction (average en¬ 
closed energy fraction is 85%). We set the sky back¬ 
ground to an annulus of a 10 arcsec width (or 13 arcsec) 
away from the center when Ag < 5 arcsec (or Ag > 5 arc¬ 
sec). We grouped all spectra with a minimum of 1 count 
per bin and applied the Cash statistic (Cash 1979) with 
the assumption of the Poisson distribution to achieve 
spectral fitting for low-count spectra. XSPEC command 
statistic cstat applies the W statistic (Wachter et al. 

1979) for the unmodeled background spectra. We per¬ 
formed spectral fitting by applying the photoelectric ab¬ 
sorption model and power law (phabs*pow in XSPEC) 
in the rest-frame energy range between 0.5 and 8.0 keV. 

Tab. 5 gives a detailed prescription for the spectral fitting 
procedure. Tab. 6 describes the X-ray properties mea¬ 
sured on these sources, and Eig. 10 shows an example 
result. 

Apart from AGNs, there are other luminous X-ray 
sources such as high-mass X-ray binaries, young super¬ 
nova remnants, and hot ionized interstellar medium re¬ 
sulting from star formation activity. In-depth studies 
have proposed a correlation between the global star for¬ 
mation rate (SER) and X-ray luminosity across the whole 
galaxy (Grimm et al. 2003; Ranalli et al. 2003; Grimes 
et al. 2005; Lehmer et al. 2010; Mineo et al. 2011). In 
particular, Ranalli et al. (2003) showed the following lin¬ 
ear relationship between 2—10 keV X-ray luminosity and 
SER based on 23 bona fide star-forming galaxies. 

SER[M 0 yr-i] = 2.0 x 10 - 4 °L 2 -iokeV (I) 



I-'X.Chandra(a-lOkeV) [erg s"‘] 

Fig. 11. — Comparison of 2—10 keV X-ray luminosity be¬ 
tween measurements (Lx Chandra) SFR-inferred estimates 

(^x, SFr)- Newly found type 1 AGNs and SpecClass 3 sources 
are shown as blue triangles and red circles, respectively. The ob¬ 
jects having more than 10% of photon pileup fraction are shown 
with filled blue triangles (3 out of 29) and filled red circles (16 out 
of 55). The black cross at the bottom right indicates the mean error 
for all sources. The central dashed line is a one-to-one fiducial line. 
Two additional dotted lines are shown to facilitate comparison. 

We can also estimate SFR by using Ha narrow emission 
line strength following Kennicutt (1998), and then de¬ 
rive the 2—10 keV X-ray luminosity (Lx, sfr) by using 
this equation. Fig. 11 compares Lx, sfr and Lx, Chandra 
for the 84 type 1 AGNs for both the newly found type 1 
AGNs (blue triangles) and the SDSS SpecClass 3 sources 
(red circles). Fluxes (horizontal axis) are calculated from 
the best-fit model. The Chandra X-ray luminosities 
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Fig. 12.— BPT diagnostic diagram for type 1 AGNs (red dots). 
Emission-line galaxies with A/N cut (> 3 for [Nil] A6584, Ha, 
[Olll] A5007, and H/3. N!^180,000) chosen from the entire OSSY 
catalogue are shown with filled contours. Dashed contour indicates 
68% distribution for type 1 AGNs. 

Chandra) are 2—4 Orders of magnitude higher than 
the SFR-inferred estimates (Lx, sfr), robustly confirm¬ 
ing that our new candidates are indeed AGNs. 

Regarding photon pileup fraction of our Chandra 
sources, we performed an investigation based on count 
rate. We explicitly marked the sources with filled sym¬ 
bols in Fig. 11 if the object has count rate per second 
greater than 0.07 which corresponds to a photon pileup 
fraction greater than 10%. Loss of photon count under¬ 
estimates X-ray luminosity (horizontal axis of Fig. 11). 
There are 19 sources (out of 84) that show photon pileup 
fraction greater than 10%, and their intrinsic X-ray lu¬ 
minosity would be found further away from the fiducial 
line in Fig. 11. Hence, photon pileup does not change our 
conclusion that 84 X-ray sources are powered by AGN. 

2.6. Type 2 AGN selection 

Based on the OSSY catalogue (N = 664,187, 2 < 0.2), 
we chose the type 2 AGNs from the Seyfert region of the 
BPT AGN diagnostic diagram (Baldwin et al. 1981). We 
used the theoretical maximum starburst model of Kew- 
ley et al. (2001) and the empirical star formation curve 
derived by Kauffmann et al. (2003) as demarcation lines. 
For the discrimination of low-ionization nuclear emission¬ 
line regions from Seyfert AGNs, we used the empirical 
demarcation of Schawinski et al. (2007). In addition, 
we used the criterion A/N > 3 as a threshold to se¬ 
cure statistical significance for the narrow emission lines 
([N il] A6584, narrow Ha, [Oiii] A5007, and H/3) used in 
the diagnostics (Fig. 12). As a result, we selected 10,108 
type 2 AGNs at 2 < 0.2. The BPT AGN diagnostics 
diagram for type 1 AGNs is presented in Fig. 12. Over¬ 
all distributions of the type 1 and type 2 AGNs versus 
redshift are shown in Fig. 13. 



redshift 


Fig. 13. — Histograms of type 1 (thick line) and type 2 AGNs 
(dark gray area) with respect to redshift; the type 1 AGN samples 
include 3,718 SpecClass 3 sources (light gray area). The hatched 
are represents the newly found type 1 AGNs. 

3. TYPE 1 AGN FRACTION 

By combining both types of AGN, we obtained the type 
1 AGN fraction, namely the number ratio between type 1 
AGNs and all AGNs, as a function of [O iii] A5007 lumi¬ 
nosity and bolometric luminosity (Fig. 14, Tab. 7). The 
type 1 AGN fraction observed (blue crosses in Fig. 14) 
increases with L[0 ill]. This trend is generally consistent 
with the prediction based on the receding torus model, 
but does not match it closely: the luminosity dependency 
of the type 1 AGN fraction is steeper at low L[Oiil] 
and becomes flattened at high L[Oiii] bins compared to 
the “standard” model (black solid curve) that was de¬ 
rived from the fit to the type 1 AGN fraction from the 
SDSS DR4. Simpson (2005) noted that the degree of 
match (to the old data) improved when the height of the 
torus was modified to depend on the nucleus luminosity 
in the form oi h oc (black dashed line) with ^ > 0. 
In the “modified” scheme, the sublimation radius of a 
dust torus gradually increases with the nuclear luminos¬ 
ity L[Oiil]. In this framework, the type 1 fraction is 
described as fi = Ni/Ntot = 1 — [1 + 3(L/Lo)^“^^]““, 
where a is set to be 0.5. In Simpson’s formulation, the 
best fit is to reproduce fi=a for a combination of (Lg, 
^); this yielded Lq = 42.37erg s“^ and ^ = 0.23 for the 
reduced xt of 1.5. The modihed model does not seem 
consistent with our new data, however. We will revisit 
this issue in the next section. 

Ghallenges remain in accurately defining the true type 
1 fraction. For example, as discussed in previous sections, 
our type 1 selection scheme probably excludes some true 
type 1 AGNs as a consequence of finding robust type 1 
AGNs. In this sense, our type 1 fraction may be con¬ 
sidered a lower limit. On the other hand, the number of 
type 2 AGNs is also uncertain. Most notably, we have se¬ 
lected type 2 AGNs only from the Seyfert category based 
on the BPT diagnostics (see Section 2.5). But, it is likely 
that there are some type 2 AGNs hidden in the composite 
and even the star-forming categories. We now attempt 
to discover them and to include them in our estimation 
of the type 1 fraction. 
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logLb„,[erg s"'] 

42.5 43.5 44.5 45.5 46.5 



logL[Oni][erg s'‘] 

Fig. 14.— Type 1 AGN fraction with respect to [Oiii] A5007 
luminosity. Bolometric luminosity is shown on the top axis, follow¬ 
ing Lbol ~ 3500L[Olll] (Heckman et al. 2004). Blue crosses and 
error bars show the observed type 1 AGN fraction. Red crosses 
and error bars show the type 1 AGN fraction estimated by sta¬ 
tistically correcting for the dilution of observed type 2 AGNs due 
to star formation. The gray band shows the range of L[Olll] of 
star-forming galaxies sampled so that their stellar mass distribu¬ 
tion matched that of type 2 AGNs. Vertical dotted lines denote the 
median and 95% distributions. The black solid curve represents the 
standard receding torus model proposed by Simpson (2005). The 
black dashed curve indicates the modified receding torus model 
which allows to vary the height of the torus according to the nu¬ 
cleus luminosity with ^=0.23. The best fits for the observed (blue) 
and dilution-corrected (red) data points are shown as blue and red 
dashed lines, including Note that the error bars are barely vis¬ 
ible due to their small magnitudes, except for the highest L[Olll] 
bin. 


3.1. Dilution correction on optical lines of type 2 
AGNs 

Our analysis starts by setting up two groups of galax¬ 
ies at z < 0.2: “emission-line” galaxies {N ^ 200, 000) 
and star-forming galaxies (TV ~ 150,000) from the OSSY 
catalog. Emission-line galaxies obviously include star¬ 
forming galaxies. We use only the galaxies with A/N 
> 3 for all four lines of [N il] A6584, Ha, [O ill] A5007 
and H/3 that are used in the BPT diagnostics. Assum¬ 
ing that the amount of emission fluxes that come from 
star formation depends mainly on its stellar mass, we 
matched every emission-line galaxy with a BPT-selected 
star-forming galaxy with the same stellar mass within a 
5% tolerance level, based on the stellar mass listed in the 
MPA-JHU catalog. The size (number of galaxies) of the 
pool of galaxies with (roughly) the same stellar mass for 
each galaxy used in our test, which depends strongly on 
stellar mass itself and the shape of stellar mass function, 
is roughly 3,400 (median). 

This match between an emission-line galaxy and a 



log([NII] X 6583/Ha) 

Fig. 15.— BPT diagnostic diagram with 5 randomly chosen ex¬ 
amples. All emission line galaxies at 2 < 0.2 having A/N > 3 on 
[N ll] A6584, Ha, [O lll] A5007, and H/9 are represented by filled con¬ 
tours. Composite objects (green stars) are shown that were moved 
to the Seyfert region (red circles) by the procedure described in the 
main text. 

star-forming galaxy solely based on stellar mass is in¬ 
valid when it is performed on a single galaxy, because 
even for the same stellar mass of galaxies, star formation 
rates can vary widely, resulting in different emission line 
strengths. However, when this exercise is applied to nu¬ 
merous galaxies, and as long as there is no serious bias 
of star formation rate in any sense in the sample, the 
method can be considered largely valid. We subtracted 
the emission line strength of the star-forming galaxy from 
the emission-line galaxies, for all four narrow lines, and 
we adjusted the position of each galaxy in the BPT di¬ 
agram based on its residual line strengths. Note that 
we have done this on all emission-line galaxies including 
star-forming galaxies themselves. 

After removing star-forming components in the emis¬ 
sion lines, 550 star-forming galaxies were moved to the 
Seyfert region of the BPT diagram, as shown in Fig. 15. 
This does not mean that there is a significant AGN com¬ 
ponent in all of these star-forming galaxies. Instead, this 
arose from the randomness of the choice of our star¬ 
forming galaxy as a subtraction. Likewise, after this 
subtraction, many star-forming galaxies ended up hav¬ 
ing negative emission lines, again as a result of the ran¬ 
domness of the process. Hence, we decided to ignore this 
result. 

The light emitted by “composite” galaxies is believed 
to include the contributions of both star-forming and 
AGN activities, and thus this exercise makes more sense 
when applied to these galaxies; applying this process 
moved 1,374 composite galaxies from the composite re¬ 
gion to the Seyfert region (Fig. 15). Again, some of 
these moves may be a result of applying an unfair star- 
formation correction due to the randomness of our pro- 
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TABLE 7 

The observed type 1 AGN fraction with the dilution corrected data 


logibol 
[erg s-l] 

logL[0 III] 
[erg s-l] 

f(type 1 AGNs/total) 
[%] 

f(type 1 AGNs/totab^-corrected 
[%] 

Ntype 1 AGN 

Ntype 2 AGN 

Tvjcorrected^ 
^^type 2 AGN 

42.5-43.5 

39.0 - 40.0 

10.90 + 0.98 

7.72 + 0.71 

110 

899 

416 

43.5 - 44.0 

40.0 - 40.5 

23.08 + 0.86 

19.60 + 0.75 

553 

1843 

426 

44.0 - 44.5 

40.5 - 41.0 

34.71 + 0.67 

32.42 + 0.63 

1778 

3344 

363 

44.5 - 45.0 

41.0 - 41.5 

42.92 + 0.71 

41.86 + 0.70 

2084 

2772 

123 

45.0-45.5 

41.5 - 42.0 

45.08 + 1.16 

44.83 + 1.15 

833 

1015 

10 

45.5 - 46.5 

42.0 - 43.0 

52.09 + 2.64 

51.80 + 2.63 

187 

172 

2 


^ Dilution corrected type 1 AGN fraction which is described in the Section 3.1 
^ Number of dilution corrected type 2 AGNs. 


cess. This yielded a total of 11,482 (10,108 + 1,374) type 
2 AGNs. The new type 1 AGN fraction that resulted 
from this “dilution” correction is shown in Fig. 14, using 
red crosses and error bars. 

The analysis is not complete, because we randomly 
chose and considered emission line fluxes of star-forming 
galaxies that only matched the stellar mass. Neverthe¬ 
less, we can expect a general trend of star-formation- 
caused dilution of AGN relying on the large-number 
statistics. The finding of new type 2 AGNs reduces the 
type 1 AGN fraction, especially in low-luminosity bins, 
but it seems clear that the current versions of the reced¬ 
ing torus model still do not reproduce the “corrected” 
type 1 fraction. 

We now attempt to find fits to our new data by adopt¬ 
ing the general form of Simpson’s formula, but treating a, 
Lq and ^ as free parameters. We assume ^ > 0 because 
this condition is essential to the receding torus model. 
The best fits were found with a = 0.12, Lq = 39.7, 
^ = 0.0, and xl = 6.6 for the “observed” data set and 
with a = 0.14, Lq = 40.1, ^ = 0.0, and xl = 8.3 for 
the “corrected” data set. The new fits are shown as blue 
and red dashed lines in Fig 14. However, these two best 
fits have goodness of fit values (xl) that are too large to 
accept, and are much larger than that found by Simpson 
for the old data (1.5). In our test, xl became smaller 
with decreasing Our best fit was found with ^ = 0. 
This implies that the introduction of “modification” (i.e., 
a positive luminosity dependence of torus height) seems 
unsupported. In fact better fits were achieved with a 
negative value of but the improvement was marginal. 
The new data call for a new and improved explanation. 

4. SUMMARY & DISCUSSION 

We have presented a process for finding new type 1 
AGNs at z < 0.2 from the SDSS DR7 database. Adding 
the new type 1 AGNs (1,835) to those in the previous 
database (3,718) that also pass our selection criteria in¬ 
creases the size of the type 1 AGN database by 49%. 
Based on the new database, we calculated a new type 1 
AGN fraction, which provides a useful test and constraint 
to the theoretical understanding of the AGN geometry. 
We are pleased to release the spectroscopic properties of 
all type 1 AGNs, given in Tab. 9, and our calculations of 
the type 1 AGN fraction based on Z/[Oiii] and the Ed¬ 
dington ratio, given in Tab. 7 and Tab. 8, respectively. 
It is our wish that these tables will be useful to the com¬ 
munity. 

We showed that our updated type 1 AGN fraction is 
not consistent with the current versions of the receding 


torus model. The type 1 AGN fraction in the low and 
intermediate luminosity bins is significantly higher than 
that given by previous data and models, even after cor¬ 
recting for dilution effects on emission lines caused by 
star formation. We found fits to the new data based on 
the basic formulation of the receding torus model, allow¬ 
ing large differences in the parameterization process, but 
these fits were poor. Thus, the new fits do not seem 
to support the hypothesis that torus depth depends on 
AGN luminosity. 

We present the type 1 AGN fraction with respect to 
black hole mass and bolometric luminosity in Fig. 16. 
We used the same bolometric correction factor, 3,500, 
to both types of AGN following Heckman et al. (2004). 
The bolometric correction factor may be uncertain; but 
if there is no serious type bias on it, our results and con¬ 
clusions in the relative sense would hold good. The color 
grids in the main panel of the Fig. 16 show that there 
is a favored area in which the type 1 AGN fraction is 
high, for higher black hole mass and bolometric luminos¬ 
ity. This ridge-shaped distribution suggests that neither 
black hole mass nor bolometric luminosity solely deter¬ 
mines the type 1 AGN fraction. Rather, the type 1 AGN 
fraction suggests a rising and falling distribution versus 
the Eddington ratio (see small panels in Fig. 16). If this 
trend is real, it may imply that the opening angle ini¬ 
tially widens with increasing Eddington ratio until the 
Eddington limit, and diminishes thereafter, which is not 
expected based on the basic scheme of the receding torus 
model. Admittedly, however, the goodness of fits (xl) is 
too poor to draw any definite conclusion yet. 

The structure of the dust torus is a key ingredient of 
understanding the AGN unification theory. However, its 
physical properties and the detailed geometry are still 
under debate. Krolik & Begelman (1988) proposed that 
the dust and gas in the torus take the form of optically 
thick, self-gravitating clouds. A clumpy and obscuring 
torus requires a small volume filling factor and a suf¬ 
ficient scale height compared to the mean freepath of 
clouds (Hdnig & Beckert 2007). In the framework of a 
clumpy torus model, the largest clouds in the torus be¬ 
come gravitationally unbound from the high-luminosity 
AGN when the AGN accretes and radiates close to the 
Eddington limit. Based on this argument, the relation¬ 
ship H/R ^ was derived where H is the scale 

height, R is inner radius of the torus, and L is the AGN 
luminosity. Note here that the scale height inversely cor¬ 
relates with the AGN luminosity unlike in the receding 
torus model. Elitzur & Shlosman (2006) described the 
origin of dust clouds as an outflow wind released from an 
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Fig. 16.— Type 1 AGN fraction with respect to bolometric luminosity and black hole mass. The main panel presents the type 1 AGN 
fraction by using a color grid to denote bins including more than 10 objects in total (type 1 AGNs + type 2 AGNs). Black dotted curves 
indicate contours (68% and 95%) for type 1 AGNs, and dashed lines indicate various Eddington ratios. (Right panels) Type 1 AGN 
fraction versus Eddington ratio for various ranges of bolometric luminosity. (Bottom panels) Type 1 AGN fraction versus Eddington ratio 
for various ranges of black hole mass. Red curves in the right and bottom panels show log-normal fits. Tab. 8 lists all data points, errors, 
and fitting results. 


TABLE 8 

Type 1 AGN fraction as a function of Eddington ratio 
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accretion disk. Because outflow velocities are expected 
to decrease with R, the torus outflow scenario also pre¬ 
dicts luminosity dependency of the type 1 AGN fraction. 
Again, however, our data shown in Fig. 16 do not seem to 
show a monotonic luminosity dependence. Rather, the 
covering factor of the torus seems to vary with the Ed¬ 
dington ratio, showing a favored spot in the black hole 
mass and bolometric luminosity diagram. 

From the perspective of black hole mass, i.e., horizon¬ 
tal axis of Fig. 16, smaller black hole mass regime shows 
lower type 1 AGN fraction. It can clearly be seen in 
Fig. 17 that type 1 AGN fraction gradually increases 
with black hole mass. As lower mass black hole allows 
more orbits close to the centre, the accretion disk is more 


heated up overionising broad Ha emitting BLRs. How¬ 
ever, it should be noted that further detailed analysis is 
required to confirm this interpretation. 

It is also noteworthy that there are discussions regard¬ 
ing the effect of the Eddington ratio and luminosity of 
AGN on the presence of BLRs. Nicastro et al. (2003) pre¬ 
sented evidence that hidden BLRs are found above the 
certain threshold value of Eddington ratio based on their 
earlier model (Nicastro 2000) and a handful of spectropo- 
larimetric samples with X-ray data. They claimed that 
BLRs are formed by accretion disk instabilities, and the 
critical radius becomes smaller than the innermost sta¬ 
ble orbit under the low enough accretion rates. Similarly, 
Elitzur & Ho (2009) showed that type 1 AGNs are not 
found below a certain luminosity based on disk-wind sce¬ 
nario. Their data supports the disappearance of BLRs at 
radiatively inefficient accretion regime which is far more 
lower than our data presented. 

In conclusion, our new type 1 AGN catalogue seems to 
pose a serious challenge to the current favorite versions 
of the AGN unified model, but the improved database 
will hopefully help us better understand the AGN clas¬ 
sifications in the end. 
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TABLE 9 

Properties of type 1 AGNs'^ 


SDSS ObjID 

RA‘> 

[deg] 

DEC‘> 

[deg] 

z 

logL[0 III] 
[erg s-l] 

EWHMbHc 
[km s“l] 

EWbHc 

[A] 

logibHa 
[erg s-l] 

log(ALBH/M0)<^ 

logLbol'^ 
[erg s-l] 

notes® 

587728309632106568 

152.29813 

2.74760 

0.06319 

39.98 

6671.6 ±89.2 

56.9 

41.18 

7 +0.05 

t -oo _o.05 

43.52 

N 

587729229300760654 
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52.48399 

0.06319 

39.85 
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18.6 

40.77 

6-83 tHs 

43.40 

N 
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-3.50442 

0.06320 

41.46 
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42.10 
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/.DZ _o 05 

45.00 

3 
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39.90373 
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Note. — (This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form 
and content.) 

^ This table is also available in the OSSY web site (http://www.gem.yonsei.ac.kr/ossy/)) 

^ J2000 coordinates 

^ Black hole mass derived from the method developed by Greene & Ho (2005) 

^ Bolometric luminosity derived from the method developed by Heckman et al. (2004) 

® Flag presenting the newly identified type 1 AGNs (‘N’) with the Chandra X-ray data (‘X’). ‘3’ indicates the SDSS SpecClass 3 objects. ‘H’ and ‘G’ indicate 
the object also identified by Hao et al. (2005) and Greene Ho (2007), respectively. 
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